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Molecular Clouds as Cosmic-ray Barometers 
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Abstract 

The advent of high sensitivity, high resolution 7-ray detectors, together with a knowledge of the dis- 
tribution of the atomic hydrogen and especially of the molecular hydrogen in the Galaxy on sub-degree 
scales creates a unique opportunity to explore the flux of cosmic rays in the Galaxy. We here present the 
new data on the distribution of the molecular hydrogen from a large region of the inner Galaxy obtained 
by the NANTEN Collaboration. We then introduce a methodology which aims to provide a test bed for 
current and future 7-ray observatories to explore the cosmic ray flux at various positions in our Galaxy. 
In particular, for a distribution of molecular clouds, as provided by the NANTEN survey, and local cosmic 
ray density as measured at the Earth, we estimate the expected GeV to TeV 7-ray signal, which can then 
be compared with observations and use to test the cosmic ray flux. 
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1. Introduction 

The question of the distribution of the cosmic rays 
(CRs) and of the CR sources in the Galaxy was addressed 
extensively by many authors (see for instance Ginzburg 
& Syrovatskii 1964 ). Secondary CR data tell us that 
CR protons and nuclei, that are accelerated by Galactic 
sources, have a confinement time in the Galaxy of about 
10 7 years at GeV energies, which decreases with increas- 
ing energy. During this time the particles accelerated by 
individual sources mix together, lose memory of their ori- 
gin, and contribute to the bulk of Galactic CRs known as 
the cosmic ray sea or CR background. 

Direct measurements provide the CR spectrum and flux 
only in the vicinity of the Solar System and the level of 
the CR spectrum and flux in other regions of the Galaxy 
is unknown. In order to investigate the CR level in dis- 
tant regions of the Galaxy and the distribution of Galactic 
CR sources one can start assuming that the CR flux and 
spectrum measured locally is representative of the typi- 
cal CR flux and spectrum present throughout the Galaxy, 
except in the proximity of isolated CR sources. The impli- 
cations of such an assumption are firstly, that there exists 
a definable CR sea of homogenised and isotropised CRs, 
far from any of the kpc-scale uniformly and continuously 
distributed CR sources, which is the typical CR flux and 
spectrum in the Galaxy; and secondly, that the Earth sits 
in a typical region far from any source where this sea can 
be measured directly. Although the 10 7 years diffusion 
time of CRs in the Galaxy seems to imply the existence of 
the sea, in principle there could indeed be a sea of CRs in 
the Galaxy, but the Earth could be located close to a so far 
undetected CR source, so that the locally measured flux is 



therefore substantially higher than the sea and the spec- 
trum substantially harder, and the local values are there- 
fore not representative of the typical case. Alternatively, 
if the CR sources are distributed quite densely, on sub- 
kpc scales in the Galaxy, there may not really be a sea CR 
spectrum at all. The local values of CR flux and spectrum 
would then vary strongly depending mostly on the prox- 
imity to the CR sources in any particular region. In this 
case the flux value at the Earth could be lower or higher, 
and the spectrum softer or harder, than any given distant 
region. Finally, it could be that the CR sources are so 
densely distributed in the Galaxy, that there might exist 
strong gradients in the CR spectrum, depending mostly 
on the proximity to the CR sources in any particular re- 
gion, and the Earth might by chance be located unusually 
far from any source and thus the locally measured flux and 
spectrum are lower and softer than typical values found 
elsewhere. 

Whilst diffusing through the Galaxy, CRs interact with 
ambient matter (both molecular and atomic) and produce 
7-rays through inelastic collisions. This emission has long 
been recognised by many authors as a unique probe of the 
parent CR distribution because the resultant 7-ray emis- 
sivity is a function only of the matter density and the 
CR spectrum (Ginzburg & Syrovatskii 1964 ; Puget & 
Stecker 1974 ; Casse & Cesarsky 1976 ; Bloemen et al 1984 
; Bloemen et al 1986 ; Lebrun et al. 1982 ; Bertsch ct al. 
1993 ; Hunter et al. 1997 ; Strong et al. 2004 ). For in- 
stance Bloemen et al 1986 investigated the galacto-centric 
distribution of cosmic rays. They defined four intervals in 
the galacto-centric radius, for which they calculated the 
column densities of atomic and molecular hydrogen. By 
comparing the 7-ray intensities with the predicted emis- 
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sivities, they derived an exponential distribution for the 
cosmic rays in the Galaxy. In the paper by Bertsch et al. 
1993 the CR flux used for calculating the 7-ray flux is var- 
ied position to position by multiplying by a "cosmic ray 
enhancement factor" c(r,l,b), which is derived by convolv- 
ing the matter distribution with a Gaussian distribution 
whose width is an adjustable parameter r . The cosmic 
ray enhancement factor c(r,l,b) takes a value ranging from 
0.1 to 1.6. The CR flux in the Galaxy varies significantly 
to that measured in the solar neighbourhood. The most 
comprehensive numerical model, which includes particle 
production and propagation in the Galaxy, is GALPROP 
(Strong et al. 2004 ) 1 . GALPROP is designed to perform 
cosmic-ray propagation calculations for nuclei, antipro- 
tons, electrons and positrons, and computes diffuse 7-rays 
and synchrotron emission. In GALPROP the CR prop- 
agation equation is solved numerically on a spatial grid, 
either in 2D with cylindrical symmetry in the Galaxy or 
in full 3D. The numerical solution proceeds in time until a 
steady-state is reached. Normalisation of protons, helium 
and electrons to experimental data is provided (all other 
isotopes are determined by the source composition and 
propagation). 7-rays and synchrotron are computed using 
interstellar gas data (for pion-decay and bremsstrahlung) 
and the interstellar radiation fields (ISRF) model (for in- 
verse Compton) . Spectra of all species on the chosen grid 
and the 7-ray and synchrotron sky-maps are output in a 
standard astronomical format for comparison with data. 
The approach adopted by Strong et al. 2004 to include 
effects of CR flux variations in the Galaxy is to calculate 
the propagation effect from CR sources, which results in 
variation of CR flux in our Galaxy (Strong & Moskalenko 
1998 ). The simulated 7-ray flux is then the background, 
which must be subtracted to recognise the sites of concen- 
trated excess 7-ray emission in a detector data, such as in 
Fermi data. 

Interestingly, based on EGRET data (Hunter et al. 
1997 ) the 7-ray emission from the inner Galaxy (300° < 
I < 60° and -10° < b < 10°, as defined in Hunter et al. 
1997 ) at energies greater than about 1 GeV exceeds, by 
some 60%, the intensity predicted by model calculations 
Bertsch et al. 1993 ; Hunter et al. 1997 . Recent ob- 
servations, performed by the LAT instrument on board 
the Fermi satellite, show that the spectra of the Galactic 
diffuse emission at Me V- GeV energies, at least at inter- 
mediate latitudes (Abdo et al. 2009 ), can be explained 
by cosmic-ray propagation models based on local obser- 
vations of cosmic-ray electron and nuclei spectra. On the 
other hand, at TeV energies the spectral features of the 
7-ray emission detected by HESS from the Galactic centre 
(GC) region (Aharonian et al. 2006 ), and the emission 
measured by Milagro from the Cygnus Region (Abdo et al. 
2007 ; Abdo et al. 2008 ), suggest that the CR flux might 
significantly vary in the different locations of the Galaxy. 

In this paper we seek to compare the CR flux and spec- 
trum observed at the Earth to the CR flux and spectrum 
in specific distant regions, derived from observations of 

1 GALPROP is available from: http://galprop.stanford.edu/ 



7-rays emitted by CRs propagating through molecular 
clouds. We first present not publicly available CO ob- 
servations obtained with the NANTEN telescope for a 
significant and interesting part of the galactic plane. We 
then demonstrate that the NANTEN data in combination 
with the Fermi data allow meaningful studies of devia- 
tions of the cosmic ray flux from the average one. The 
level of the sea of lower energy CRs (< few GeV) has 
been traced by the ionisation rate inferred from differ- 
ent molecular species seen in molecular clouds (Indriolo 
et al. 2007 ; Indriolo et al. 2009 ). The level of the CR 
sea above 1 GeV will then be inferred by combining the 
7-ray emission from molecular clouds with the data from 
the Leiden/ Argentine/Bonn (LAB) survey of the Galactic 
atomic hydrogen (Karberla et al. 2005 ) and from the 
NANTEN survey of Galactic molecular hydrogen (Fukui 
et al 1999 ; Mizuno & Fukui 2004 ; Fukui et al. 2001 
; Fukui et al. 2008 ). Studying the 7-ray emission from 
single molecular clouds the differential, local distribution 
of the cosmic ray flux can be inferred. The detection of 
even a single under-luminous cloud with respect to the 
predictions based on the local CR density would cast se- 
vere doubt on the assumption that the level of the ob- 
served CR flux is representative of the average CR flux 
and spectrum in the Galaxy. 

In Section 2 we will present new data on the distribution 
of the molecular hydrogen from the region of the inner 
Galaxy, which spans Galactic longitude 340° < I < 350° 
and Galactic latitude —5° < b < 5°, obtained by the 
NANTEN Collaboration. We will also briefly review the 
results of the LAB survey of the Galactic atomic hydro- 
gen. A description of our approach will be given in Section 

3. The method will be here applied to the region of the 
Galactic Plane described above to provide a 'proof-of- 
concept'. A more comprehensive study on a larger scale 
will be discussed in a future paper. Some predictions con- 
cerning the 7-ray and cosmic ray fluxes from molecular 
clouds including their uncertainties are given in Section 

4, where we will also discuss the observational prospects 
for present and future observatories. Our conclusions are 
given in Section 5. 

2. The gas distribution 

2.1. The NANTEN survey of molecular hydrogen 

The NANTEN instrument is a 4 m millimetre/sub- 
millimetre telescope which has surveyed the southern sky 
using the 12 CO (J=l-0) emission line at 115.271 GHz 
(A = 2.6 mm). The angular resolution at this frequency is 
4', with a mass sensitivity of about 100 M Q at the Galactic 
centre (GC), assuming a distance of 8.5 kpc for the GC. 
The NANTEN data set has 1 km/s resolution in velocity 
(Mizuno & Fukui 2004 ; Fukui et al 1999 ; Fukui et al. 
2001 ; Fukui et al. 2008 ). The conversion factor, X, 
which is used to translate the intensity of CO molecular 
line emission into H 2 column density (Nakanishi & Sofue 
2003 ; Nakanishi & Sofue 2006 ), is assumed to be 

V = 1.4 x 10 20 e Wllkpc) [cm^K^km^s], (1) 
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Fig. 1. The radial velocity vs heliocentric distance estimated 
from flat rotation of the galactic disk is shown. The solid line 
is at (l,b) (340,0) and the dashed line is at (l,b) (350,0). 

where R is distance to the GC. This conversion fac- 
tor, which is based on the metalicity distribution in the 
Galaxy, is a function of the distance to the GC. The un- 
certainty in the determination of the distance to the GC 
and in the determination of the calibration itself leads to 
an uncertainty in X of a factor of 50% (Arimoto et al. 
1996 ). 

The final data product for the region, which spans 
Galactic longitude 340° < I < 350° and Galactic latitude 
—5° < b < 5°, consists of a 150 x 150 x 600 element three- 
dimensional image of the molecular hydrogen distribution 
as a function of Galactic longitude, latitude and veloc- 
ity. Assuming a flat rotation curve model of the Galaxy 
with uniform velocity equal to 220km/s, the data-cubes 
in longitude, latitude and velocity are transformed to 
data-cubes in longitude, latitude and heliocentric distance 
(Nakanishi & Sofue 2003 ; Nakanishi & Sofue 2006 ). 

Figure 1 shows a plot of the radial velocity versus helio- 
centric distance used to convert the velocity information 
into distance. There are two distance solutions for each 
radial velocity in the inner Galaxy. In order to resolve this 
near-far distance ambiguity, the method by (Nakanishi & 
Sofue 2003 ; Nakanishi & Sofue 2006 ) is used. The prin- 
ciple uncertainty in the determination of the distances in 
the NANTEN data-set comes from errors in the accuracy 
of the radial velocity estimates. The combination of tur- 
bulence and peculiar proper motions of the gas causes 
spectral line broadening, which introduces uncertainties 
in the radial velocity, that can be as large as 15 km/s (a 
giant molecular cloud typically shows velocity width of 10 
km/s, and the cloud-cloud motion is 3-3.9 km/s (Alvarez 
et al. 1990 ; Clemens 1985 )). This corresponds to an 
uncertainty of less than 2 kpc in the distance. The de- 
pendence of the uncertainties in the radial velocity on the 
particular location is small for our dataset, since the lon- 
gitude range used here includes disk gas, and excludes the 
Galactic centre and molecular ring. For disk gas, it is, in 
fact, unlikely to have significant systematic variation in ei- 
ther cloud-cloud velocity dispersion or typical line widths 



with location. Correspondingly, the distance error arising 
from molecular gas velocity dispersion is likely to be fairly 
constant throughout our data sample. 
2.1.1. Gas morphology 

Figure 2 (Left) and 2 (Right) show NANTEN observa- 
tions of the 12 CO (J=l-0) emission line from a large area 
of the Galactic plane. This segment spans Galactic longi- 
tude 340° < I < 350° and Galactic latitude -5° < b < 5°. 
Figure 2 (Left) shows the integrated emissivity in units 
of cm -3 km s" 1 from this region. Inspection of Figure 
2 (Left) shows that the molecular material is predomi- 
nantly confined to < 2° of the plane. The longitude span 
places this region in between the Galactic centre and the 
Norma arm (Duncan et al. 1995 ) so that there should be 
no contamination from depth effects from looking down an 
arm. This is reflected in the distribution of heliocentric 
distance in the region in Figure 2 (Right), which shows 
the distance as a function of position. Figure 2 (Right) 
shows, in fact, a moment 1 map of the region, which is 
^[/d]/^[I], where I is the intensity and d is the dis- 
tance of the gas and has units of kilo-parsec. There are 
a few regions of emission which are very far away (cor- 
responding to the dark pixels). However, these tend to 
be isolated, small segments on the edge of the plane, and 
hence they should contribute a negligible amount to the 
7-ray emissivity. Specifically, Figure 2 (Right) shows that 
the bulk of the observed emission originates at distances 
of < 2 kpc, whilst there is some emission, mainly confined 
to the central regions of the Galactic plane, which are at 
greater distances. 

2.2. The Leiden/Argentine/Bonn Galactic HI Survey 

We use data from the Leiden/ Argentine/Bonn (LAB) 
Galactic HI Survey 2 (Karberla et al. 2005 ) for the dis- 
tribution of the atomic hydrogen. The LAB observations 
were centred on the 1420 MHz (A21 cm) line with a band- 
width of 5 MHz. The velocity axis spans —450 km s _1 to 
400 km s _1 giving a final velocity resolution of 1.3 km s _1 . 
The data combines observations from three telescopes at 
a common resolution of 0.6° with a final sensitivity of 
~ 0.1 K. We used the same procedure described in 2.1 to 
transform the velocity axis into distance. 

3. Testing the level of the CR background 

The emissivity of molecular clouds located far away 
from CR sources, so called passive molecular clouds, i.e. 
clouds which arc illuminated by the supposedly existing 
CR background, can be used to probe the level of the 
CR sea (Issa and Wolfendale 1981 ). Given that the 7- 
ray -emission from the molecular cloud depends only upon 
the total mass of the cloud, M, and its distance from the 
Earth, d, the CR flux, &cr, m the cloud is uniquely de- 
termined as 



available from: 

http://www.astro.uni-bonn.de/~webaiub/english/tools 
labsurvey . php 
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where F 7 is the integral 7-ray flux from the cloud and M = 
nV , with n the gas number density and V the volume of 
the cloud. Under the assumption that the CR background 
is equal to the locally observed CR flux, the calculated 7- 
ray flux from the cloud can be compared to the observed 
7-ray flux in order to probe the CR spectrum in distant 
regions of the Galaxy. The detection of under- luminous 
clouds with the respect to predictions based on the CR 
flux at Earth would suggest that the local CR density 
is enhanced with respect to the Galactic average density. 
This would cast doubts on the assumption that the local 
CRs are produced only by distant sources, and that the 
CR flux and spectrum measured locally is representative 
of the typical CR flux and spectrum present throughout 
the Galaxy (Aharonian 200f ). 

3.1. Diffusion of CRs into molecular clouds 

The diffusion of CRs into molecular clouds depends 
upon the highly uncertain diffusion coefficient. CRs can 
penetrate clouds if the diffusion coefficient inside the 
molecular clouds is the same as the average diffusion co- 
efficient in the Galaxy, as derived from spallation mea- 
surements. CRs are effectively excluded if the diffusion 
coefficient is suppressed (Gabici, Aharonian & Blasi 2007 
). However, it has been shown that CRs at TeV en- 
ergies can diffuse into even the densest parts of molec- 
ular clouds. Whilst GeV energy CRs can diffuse into 
typical Galactic molecular clouds (Gabici et al. 2009 ), 
they have trouble penetrating the densest parts of par- 
ticular molecular clouds, such as the GC giant molecular 
cloud, Sagittarius B2 ((Jones et al. 2009 ) and references 
therein). 

3.2. Calculational procedure 

Under the assumption that the CR background in the 
region of longitude range 340° < I < 350° and latitude 
range —5° < b < 5° is uniform and equal to the locally ob- 
served CR flux we calculate the spectral features and the 
longitudinal and latitudinal profiles of the 7-ray emission. 
We will then proceed to infer the missing information con- 
cerning the distance from the Sun to complete the 3D view 
of the 7-ray emission. 
3.2.1. Spectral feature of the emission 

The differential photon flux at Earth produced by CRs 
interacting with the gas from a given region in the Galaxy 
is 



d/V-v 



dAdE^dtdfl 



din 



do 



dE n 



n(l,b,l d )J(E p )dE p (3) 



where the integrals are over the heliocentric distance, Id 
and over the energy of the protons, E p . n(l,b,ld) is the 
gas density as function of the heliocentric coordinates, lat- 
itude, b, longitude, I, and Id, and J(E p ) is the proton 
spectrum. In Equation (3) the spectral dependence of the 
photons emitted by CR protons, expressed in terms of 
the differential cross section, do- p _ >7 /dE p , is calculated 
using the parametrisation of Kelner et al. 2006 . A mul- 
tiplication factor of 1.5, applied to the proton spectrum, 
accounts for the contribution to the emission from heavier 



nuclei both in CRs and in the interstellar medium (Dermer 
1986 ; Mori 1997 ). From Equation (3) for a uniform dis- 
tribution of CRs the 7-ray spectrum is given as 

dN~, 



dAdE^dtdfl 



J dEp ^dE^ J{Ep) J dldn ^ ld ^ W 

where the integral over the line of sight, l d , defines the 
gas column density for a given direction in latitude and 
longitude. The average atomic and molecular hydrogen 
column densities for the region of the Galaxy of longitude 
range 340° < / < 350° and latitude range -5° < b < 5° are 
3.3 x 10 21 enr 2 and 1.4 x 10 22 cur 2 , respectively The CR 
spectrum is assumed to be equal to the CR flux measured 
at the top of the Earth's atmosphere. Below 10 6 GeV we 
here adopt the value of 



J(E p ) = 1.8^ p " 2 - 7 cm- 2 s- 1 sr- 1 GeV- 



(5) 



taken from the Particle Data Group 2008 . A change in the 
spectral index of the measured CR flux is observed above 
10 6 GeV, which we do not consider as we are interested 
in the 7-ray emission below 100 TeV. 

3.2.2. Longitudinal and latitudinal profiles of the emis- 
sion 

The latitudinal and longitudinal profiles of the 7- 
ray emission due to protons scattering off the atomic 
and molecular hydrogen from the region which spans 
Galactic longitude 340° < I < 350° and Galactic latitude 
— 5°<fo<5° are shown in Figure 3 (Left) and (Right) re- 
spectively. The 7-ray profiles which the Fermi telescope 
would observe are plotted for 1 GeV and 10 GeV. For a 
Cherenkov telescope, such as CTA, we show the 7-ray pro- 
files at 100 GeV and 1 TeV. The instrument point spread 
functions at the different energies are also taken into ac- 
count. 

A peak in the emission at longitude of about 345.7° 
close to the Galactic Plane is clearly visible, next to a 
dip in the longitude profile. While the atomic gas is gen- 
erally broadly distributed along the Galactic Plane, the 
molecular hydrogen is less uniformly distributed and the 
peaks in the 7-ray longitude profiles correspond to the 
locations of highest molecular gas column density. The 
peaks in the longitudinal profile reveal the directions in 
the Galaxy where massive clouds associated with spiral 
arms are aligned along the line of sight. In the following 
discussion we will further analyse the region of latitude 
range 1° < b < 1.7° and longitude range 345.3° < I < 346°, 
where the peak in the 7-ray profiles is located. 

3.2.3. Localisation of the emission along the line of sight 
distance 

In order to localise the region along the line of sight, Id, 
where most 7-rays are emitted, Figures 4 and 5 (Middle) 
show the atomic and molecular hydrogen densities as a 
function of the line of sight distance, averaged over the 
inner 1° < b < 1.7° of the Galactic plane in the region of 
longitudes 345.3° < I < 346°. The value of the gas density 
is given every 50 pc from 50 pc up to 30 kpc from the 
Sun. The density of the atomic hydrogen as function of 
the heliocentric distance, l d , shown in Figure 4 (Middle), 
tends to be quite uniform. On the other hand, the plot 
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Fig. 2. (Left): Integrated intensity image for the CO (J=l-0) line emission obtained by the NANTEN telescope for the region 
340° < I < 350° and —5° < b < 5°. The intensity scale on the right is in units of km s _1 cm -3 with the 4' beam located in the lower 
left-hand corner. (Right): Position-distance image (i.e., moment 1) of the region of interest. The intensity scale on the right is in 
units of kpc. A logarithmic scaling has been applied to both images to emphasise the low-level emission. 
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Fig. 3. (Left 4 panels) The latitudinal profile of the 7-ray emission due to CRs scattering off atomic and molecular hydrogen 
at different energies which would be measured by Fermi and a Cherenkov telescope such as CTA from the region —5° < b < 5°, 
integrated over the longitude range 340° < / < 350°. The emission is convolved with the energy dependent PSF of Fermi at 1 GeV 
and 10 GeV, and of a future Cherenkov telescope at 100 GeV and 1 TeV. The energy dependence of the Fermi PSF can be found 
at:http://www-glast. slac.stanford.edu/software/IS/glastlatperformance. htm. We adopt a point spread function of 0.05° at 100 GeV 
and 0.02° at 1 TeV. The dotted green lines are for the emission arising from CRs scattering off molecular hydrogen, the dashed red 
ones for the emission arising from CRs scattering off atomic hydrogen and the solid black ones for the sum. The 7-ray profile is 
in units of cm~ 2 s~ 1 sr _1 TeV — 1 . (Right 4 panels) The longitudinal profile of the 7-ray emission from the region 340° < I < 350°, 
integrated over the latitude range —5° < b < 5°. The emission is convolved with the energy dependent PSF of Fermi at 1 GeV and 
10 GeV, and of a future Cherenkov telescope such as CTA at 100 GeV and 1 TeV. The dotted green lines are for the emission arising 
from CRs scattering off molecular hydrogen, the dashed red ones for the emission arising from CRs scattering off atomic hydrogen 
and the solid black ones for the sum. The 7-ray profile is in units of cm - 2 s — 1 sr _1 TeV — 1 . 
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in Figure 5 (Middle) shows that much of the molecular 
material is located within < 2 — 3 kpc, which reinforces 
the evidence presented in Figure 2 (Right). As Figure 
5 (Middle) shows, the molecular hydrogen density peaks 
to riH 2 = 20 cm~ 3 in this region. 

The total mass (H 2 and HI) in the element of volume, 
AV = ljcos(b)AbAlAl d , with A6 = AZ = 0.01 rad (=0.7°) 
and Aid = 50 pc around the directions of longitude 345.7° 
and latitude 1.3°, is 



l<b<1.7, 345.3<1<346 



AM = n(l,b,l d )AV . 



(6) 



AM is calculated every 50 parsecs along the line of the 
sight distance, l d , corresponding to each value of the gas 
density n(l,b,l d ). In Figures 4 and 5 (Top) we plot the 
mass per unit distance AM / Al d in a volume clement in 
units of M 5 Icpc^ 1 where M 5 = 10 5 M Q . 

Assuming the local CR flux, the unit volume emis- 
sivity above a given energy, E* , is proportional to the 
total mass in the element of volume, AV, and in pho- 



tons cm 



kpc 1 it is 



AF(E > E*) = 1.2 x 10 



Id 



-13 



1 kpc 



E* \ 
1 TeVy 

'AM 



-1.7 



Ms 



(7) 

The cmissivity per unit volume as a function of the line 
of sight distance, AF(E > E*)/Al d , is plotted in Figures 
4 and 5 (Bottom). Figures 4 and 5 (Bottom) show that 
most of the 7-ray emission from the region of longitude 
345.3° < I < 346° and latitude 1° < b < 1.7° is produced 
between 0.5 and 3 kpc. 

The combined longitudinal, latitudinal and distance 
profiles provide together a tomographical view of the 7- 
ray emission. The investigation of the distribution of the 
7-ray emission gives the locations in the Galaxy where 
most 7-ray emission is produced within a reasonably lim- 
ited region, where a single or at most few massive clouds 
are located, making it possible to constrain the CR flux 
in this particular region. 

3.3. Leptonic emission and secondary electrons 

3.3.1. Leptonic Emission: The Inverse Compton and 
Bremsstrahlung Contribution 
Inverse Compton scattering of cosmic ray electrons re- 
sults in a non-negligible contribution to the overall diffuse 
7-ray emission at GeV, and perhaps also at TeV ener- 
gies (see e.g. (Aharonian & Atoyan 2000 )). However the 
relative contribution of this component (compared to the 
bremsstrahlung and 7r°-decay 7-rays) from specific dense 
regions is significantly reduced because of the enhanced 
gas density. The ratio of the bremsstrahlung-to 7r°-decay 
7 rays does not depend on the density. However, the con- 
tribution of the bremsstrahlung to the total diffuse emis- 
sion above 1 GeV is rather low, so it can be ignored. This 
is true also for the bremsstrahlung from secondary elec- 
trons. Therefore in this paper we will focus only on the 
7r°-decay component. 
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Fig. 4. (Top): Total HI mass per unit distance ((AM / Ad)) 
in a volume element in 10 5 Mq kpc -1 from the direction 
345.7° in Galactic longitude and 1.3° in Galactic latitude is 
shown. (Middle): Gas density as a function of heliocentric 
distance l^. (Bottom): Unit volume emissivity above 1 GeV 
as a function of Id in units of cm - 2 s -1 kpc — . 



3.3.2. Possible Role of Secondary Electrons and 
Positrons 

When CRs interact with ambient matter through in- 
elastic collisions they produce not only neutral pions, but 
also charged pions. These charged pions emit secondary 
electrons, positrons and neutrinos. Thus, in addition to 
the 7-ray emissivity due to neutral pion decay, there will 
be a simultaneous radio synchrotron component due to 
the secondary electrons and positrons produced concomi- 
tantly with the neutral pions. Whilst the emissivity of 
such electrons is beyond the scope of this paper, we pause 
only to note that their emissivity is proportional to the 
mass of ambient matter and the CR flux in most regions. 
In the densest regions (i.e. > 10 4 cm~ 3 ) however, the 
synchrotron emissivity will scale as the volume because of 
the increased efficiency of bremsstrahlung cooling (Jones 
et al. 2009 ) . Since synchrotron emission from secondaries 
is not readily observable in the cases so far attempted 
(Jones et al. 2008 ; Protheroe et al. 2008 ), this method 
of checking the validity of the CR- sea assumptions at GeV 
energies is probably not feasible. 
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Fig. 5. (Top): Total H 2 mass per unit distance ((AM/ Ad)) 
in a volume element in 10 s Mq kpc -1 from the direction 
345.7° in Galactic longitude and 1.3° in Galactic latitude is 
shown. (Middle): Gas density as a function of distance d. 
(Bottom): Flux above 1 GeV in units of cm - 2 s — 1 kpc — 1 . 



4. Results and Discussion 

Table 1 shows the total hydrogen mass from the LAB 
and NANTEN surveys and the integral 7-ray flux above 
1 GeV and above 100 GeV as a function of the line of sight, 
l d , from the region 345°. 3 < I < 346° and 1° < b < 1° .7 
centred on I = 345°. 7 and b — 1°.3. The values for the 
masses and fluxes are obtained by integrating Equation 
(6) and (7), respectively, along given portions of the dis- 
tance, Id- It is instructive to note that 2% of the emission 
comes from the region d < 0.5 kpc, 85% from the region 
0.5 kpc < d < 3 kpc, and 13% from the region d > 3 kpc. 
The peak in the 7-ray emission in Figure 3 from the direc- 
tion 345.7° close to the Galactic plane is mostly produced 
within 0.5 kpc and 3 kpc distance from the Sun. Thus the 
7-ray emission from this direction provides a unique probe 
of the CR spectrum in 0.5-3 kpc, provided the 7 emission 
morphology is shown to match the gas in this distance 
range. In Figure 6 we show the predicted 7-ray spectra for 
the same region, 345.3° < I < 346° and 1° < b < 1.7°, aris- 
ing from the CR spectrum in Equation (5). As shown in 
Figure 7, the contribution to the emission from the atomic 
hydrogen in the 0.5kpc < d < 3kpc amounts to about 10% 
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Fig. 6. The spectra of the 7-ray emission from the region 
345.3° < I < 346° and 1° < b < 1.7° arising from the observed 
CR spectrum as in Equation (5). The atomic hydrogen col- 
umn density is 2.1 X 10 22 cm -2 and 8.0 X 10 22 cm -2 for the 
molecular hydrogen. 

of the total contribution from atomic and molecular hy- 
drogen. The emission beyond 3 kpc arises predominantly 
(80%) from the more uniformly distributed HI gas. 

The Fermi Collaboration quotes an integral sensitivity 
of 5 x 10~ 10 photons cm~ 2 s^ 1 above 1 GeV from a re- 
gion which has an angular resolution of 0.7 by 0.7 de- 
grees. This makes the detection of the predicted emission 
for Fermi possible. Future Cherenkov telescopes, such as 
CTA, will need to reach a sensitivity of about 5 xl0~ 12 
photons cm~ 2 s _1 for observing a region of 0.7 by 0.7 
degrees above 100 GeV. 

We have here presented a methodology to test the cos- 
mic ray flux in discrete distant regions of the Galaxy by 
combining the most recent high resolution 7-ray and gas 
data. By assuming a flat rotation curve model of the 
Galaxy we have estimated the gas number density at dis- 
crete positions as a function of the heliocentric latitude, 
longitude and distance. Through the knowledge of the 
gas number density we have determined the total mass 
of the gas in discrete regions of the Galaxy and there we 
have predicted the 7-ray emissivity under the assumption 
that the CR flux is equal to that measured at Earth. By 
comparing the predicted and the measured 7-ray flux we 
have obtained a test of the level of the CR flux in these 
regions. The CR flux has been calculated as a function of 
the Galactic longitude, latitude and heliocentric distance, 
providing a sort of tomography of the CR flux in different 
locations of the Galaxy. 

For our investigations we have used the higher reso- 
lution 12 CO (viz. 4' NANTEN data) and LAB HI sur- 
veys. The resolution of the LAB HI data used in the 



8 Astronomical Society of Japan [Vol. , 

Table 1. Summary of parameters as a function of distance for the region 345.3° < I < 346° and 1° < b < 1.7° 



Region 

(kpc) 

(1) 


Mass 
(Mb ) 
(2) 


Mass (H 2 ) 

(M 5 ) 

(3) 


Mass (HI) 

(M 5 ) 

(4) 


F( E > 1 GeV) 
(cm~ 2 s _1 ) 
(5) 


F(E > 100 GcV) 
(cm~ 2 s _1 ) 
(6) 


0.5 < d < 30 


20.2 


4.2 


16.0 


1.7 xl0~ s 


6.8 xlO- 12 


d < 0.5 


0.0015 


0.0005 


0.001 


2.9 xl0~ 10 


1.2 xl0~ 13 


0.5 < d < 3 


3.2 


2.8 


0.4 


1.4 xl0~ 8 


5.7 xl0~ 12 


d > 3 


17.0 


1.4 


15.6 


2.3 xl0~ 9 


9.1 xl0~ 13 




Fig. 7. Flux above 1 GeV as function of the line of sight dis- 
tance in units of cm -2 s _1 kpc -1 , produced by CRs interact- 
ing with the atomic and molecular gas. The atomic hydrogen 
is uniformly distributed in the Galactic Plane, whereas the 
molecular hydrogen is more localised. 

present investigation is 0.6 deg. The high resolution of 
the NANTEN data allows for the first time detailed com- 
parisons of our 7-ray predictions with new TeV and multi- 
GeV data. In future we will make use of the higher resolu- 
tion (2 arcmin) Southern Galactic Plane HI Survey avail- 
able at 3 , in order for our predictions to match with higher 
resolution TeV data in cases where the 7-ray component 
from HI is similar to or dominates over the 7-ray com- 
ponent from H2 (e.g from Figure 4 in the 10 to 15 kpc 
range) . 

4-1- Uncertainties in the CR flux estimate 

Uncertainties in the distance estimates affect the cal- 
culation of the CR flux in Equation (2) in three ways. 
Firstly, in the determination of the conversion factor X in 
Equation (1), by which the emissivity of the CO line is 
converted to ambient matter density. As previously dis- 
cussed, the uncertainty in the X factor is 50% (Arimoto 
et al. 1996 ). Secondly, there are uncertainties in the de- 

3 http : //www. atnf . csiro . au/research/HI/sgps/queryForm. 
html 



termination of the number density of H2 molecules from 
the measured CO intensity. Specifically, the error on the 
distance (especially at small distances; c.f. Fig.l) is as 
large as 2 kpc. We note however, that for extended sources 
(i.e., sources whose intrinsic size is larger than the beam 
size), Equation (2) shows that the CR flux density is pro- 
portional to the distance squared and inversely propor- 
tional to the mass. Since the mass scales as the distance 
square, the relative errors exactly cancel. For unresolved 
sources (i.e., sources whose intrinsic size is smaller than 
the beam size), we note that the above errors still hold, 
however, there will be a systematic underestimation of the 
true CR flux density because of beam dilution. Thus esti- 
mates for the CR flux density need to be estimated using 
extended sources only. 

It is instructive to note that only 85% of the emission, 
which would be measured from the direction 345.3° < 
I < 346° and 1° < b < 1.7°, is produced in the region 
0.5kpc < d < 3kpc. Therefore only this 85 % of the emis- 
sion is to enter in Equation (2) in order to estimate the 
CR flux in the region 0.5 kpc < d < 3 kpc. Moreover, part 
of the 7-ray emission, which would be measured from the 
region 345.3° < I < 346° and 1° < b < 1.7°, is produced 
through competing IC and non-thermal bremsstrahlung 
mechanisms of primary and secondary electrons. In par- 
ticular, in the energy range between 10 MeV and 100 MeV 
bremsstrahlung processes of background and secondary 
electrons produce a non negligible contribution to the 7- 
ray emission, so that the presented method can be applied 
to the analysis of the 7-ray emission above 100 MeV-1 
GeV, where moreover the angular resolution of the Fermi 
LAT detector is substantially better. Finally unknown 
CR source (s) might produce an enhancement of the local 
CR background. 

5. Conclusions 

The CR level in distant regions of the Galaxy and the 
distribution of Galactic CR sources are unknown. 7- 
rays emitted when CR protons and nuclei interact with 
the atomic and molecular gas have been long recognised as 
a unique probe of the CR flux and spectrum in the Galaxy 
(Ginzburg & Syrovatskii 1964 ; Puget & Stecker 1974 
; Casse & Cesarsky 1976 ; Bloemen et al 1984 ; Bloemen 
et al 1986 ; Lebrun et al. 1982 ; Bertsch et al. 1993 
; Strong et al. 2004 ). 

We have here introduced a methodology which aims to 
provide a test bed for current and future 7-ray observato- 
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ries to explore the cosmic ray flux at various positions in 
our Galaxy. In order to investigate the CR spectrum in 
distant regions of the Galaxy we have here assumed that 
the CR flux and spectrum measured at the Earth is rep- 
resentative of the typical CR flux and spectrum present 
throughout the Galaxy. The 7-ray emission arising from 
such CR spectrum has been calculated in specific regions 
of the Galaxy, using the knowledge on the distribution of 
the atomic and molecular gas, as given by the LAB and 
Nanten surveys, respectively. The peaks in the longitu- 
dinal and latitudinal distributions of the gas correspond 
to directions in the Galaxy, where only one or at most 
few massive clouds contribute to the total mass. This can 
be demonstrated by studying how the gas is distributed 
along the line of sight distance. We have then exploited 
this knowledge of the gas number density point by point 
to determine the total mass of the gas of some MCs of the 
Galaxy. From the total mass of the atomic and molecular 
gas we have calculated the level of the 7-ray emission in 
such gas clouds. The predicted 7-ray emission from such 
gas clouds can be compared with observations of present 
and future 7-ray detectors and provide a test of the CR 
density in distant locations of the Galaxy. We have, in 
particular, demonstrated that the new high quality data 
of NANTEN and the performance of the Fermi LAT tele- 
scope are sufficient for meaningful probes of CR densities 
in specific localised regions. 

Equation (2) shows that the CR flux density is propor- 
tional to the distance squared and to the 7-ray ray flux 
and inversely proportional to the mass of the MC. The er- 
ror in the determination of distances (especially at small 
distances; c.f. Fig.l) is as large as 2 kpc. We note however, 
that for extended sources (i.e., sources whose intrinsic size 
is larger than the beam size) the mass scales as the dis- 
tance square and the relative errors on mass and distance 
exactly cancel. Moreover, above 100 MeV - 1 GeV, the 
competing leptonic processes, bremsstrahlung and IC of 
background and secondary electrons, produce a negligible 
contribution to the 7-ray emission. The emission from 
molecular clouds, located far from strong CR sources, is 
then the footprint of the CR sea within the cloud and 
molecular clouds serve then as barometers of the CR pres- 
sure. The detection of even a single under-luminous cloud 
with respect to the predictions based on the local CR den- 
sity would cast severe doubt on the assumption that the 
level of the observed CR flux is representative of the aver- 
age CR flux and spectrum in the Galaxy. This conclusion 
is even reinforced if additional gas, not yet revealed, is 
present in the cloud. In this regard, by comparing various 
interstellar gas tracers, EGRET 7-ray sky maps, HI and 
CO lines but also dust thermal emission, Grenier et al. 
2005 unveiled clouds of dark gas surrounding the nearby 
CO clouds. Grenier et al. 2005 concluded that the addi- 
tional mass from dark clouds implies that the cosmic ray 
densities are overestimated in regions rich in dark clouds. 

A detailed analysis of the spectral features of the 7- 
ray emission from various regions of the Galaxy is also 
crucial. In fact, deviations from the standard spectrum 
are indications of local CR sources. The regions in the 



Galaxy where the spectral features and/or the normalisa- 
tion of the 7-ray spectrum cannot be explained by inter- 
actions of a CR spectrum similar to that measured close 
to the Sun, are of interest to probe the sites of Galactic 
CR injection, such as SNRs and pulsars. 

The Fermi observatory, which detects 7-rays in the en- 
ergy range between tens of MeVs and hundreds of GeVs 
and which benefits by a large effective area and good an- 
gular resolution, is an ideal instrument to perform sur- 
vey of extended diffuse sources such as molecular clouds 
in the Galaxy. Additionally, future VHE detectors, such 
as CTA, will offer improved sensitivity and angular res- 
olution, close to the angular resolution of the NANTEN 
survey, and could investigate the morphology of the 7- 
ray emission and compare it to the structures of molec- 
ular hydrogen clouds in the Galaxy. Future IACT tele- 
scopes ability to resolve clouds will of course strongly de- 
pend upon the instrument sensitivity to detect extended 
sources. On the other hand, future large field of view ex- 
periments, such as HAWC, will be able to investigate the 
level of the cosmic ray background on large portions of 
the sky (Aharonian, Buckley, Kifune & Sinnis 2008 ). 

Finally, we note that we have applied our analysis of the 
spectral line observations to predict the 7-ray emission in 
only one distant region. It serves mostly as a 'proof-of- 
concept', to be followed by a more comprehensive study 
of the entire Galaxy. 
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